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A flexible family of distributions on the cylinder 

Shonosuke Sugasawa*, Kunio Shimizu* and Shogo Kato* 


Abstract 

We propose a flexible family of distributions, generalized f-distributions, on the cylinder which 
is obtained as a conditional distribution of a trivariate t distribution. The new distribution has 
unimodality or bimodality, symmetry or asymmetry, depending on the values of parameters and 
flexibly fits the cylindrical data. The circular marginal of this distribution is distributed as a 
generalized t-distribution on the circle. Some other properties are also investigated. The proposed 
distribution is applied to the real cylindrical data. 

Key words and phrases: circular-linear correlation, circular-linear regression, generalized von 
Mises distribution, Johnson-We hr ly model, Mardia-Sutton model 


1 Introduction 

Directional or circular data often appear in a variety of scientific fields and various stochastic models 
have been proposed for analyzing such data. For univariate circular data, there have been many 
distributions investigated in terms of both tractability and applicability, see Jones and Pewsey (2005), 
Kato and Jones (2010) and Kato and Jones (2015). However, we sometimes encounter situations 
which involve both circular and linear variables, namely cylindrical data such as the pair of wind 
direction and temperature (Mardia and Sutton, 1978) or directions and distances of animal movements 
(Fisher, 1993). For such data, the distribution on the cylinder is needed, but there are not so many 
distributions compared to univariate circular distributions. We give a brief review below for several 
cylindrical distributions known in the literature. Johnson and Wehrly (1978) gave a distribution 
based on the principle of maximum entropy subject to constraints on certain moments, and Mardia 
and Sutton (1978) provided another distribution as a conditional distribution of a trivariate normal 
distribution or a maximum entropy distribution. An extension of the distribution by Mardia and 
Sutton (1978) was studied by Kato and Shimizu (2008), which can also be derived as a maximum 
entropy distribution or a conditional of a trivariate normal. 

In this paper, we propose the generalized t-distribution on the cylinder, which is a natural extension 
of the member of the exponential family given by Kato and Shimizu (2008). The proposed distribution 
is also regarded as a cylindrical extension of the generalized t-distribution on the circle proposed by 
Siew et al. (2008). In fact, the circular marginal distribution is the generalized t-distribution on 
the circle. The proposed distribution can be obtained as a conditional distribution of a trivariate 
t distribution and characterized as the maximum /3-entropy distribution. This is a quite flexible 
distribution which allows for both asymmetry and variations in tail weight in terms of parameter 
values. We investigate some properties such as marginal and conditional distributions, modality, 
moments, circular-linear correlation and skewness. We briefly discuss a circular-linear regression 
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model derived from the conditional distribution. For practical use, we provide an iterative algorithm 
for parameter estimation. 

Subsequent sections are organized as follows. Section 2 provides a derivation of the distribution. 
In Section 3 some properties of the new distribution are studied. We apply the proposed distribution 
to the data set of the wind direction and ozone level given in Johnson and Wehrly (1977) in Section 4. 


2 Derivation 

Suppose that a trivariate random vector W is distributed as a trivariate t distribution with degrees of 
freedom a, mean vector r/ = (771 , 772 ? ^73 ) x € M 3 and variance-covariance matrix E, where 


a 


E = 


1 

P 12 O 1 V 2 


P 12&1&2 P 13 C 1 C 3 


(To 


P23<T2<T3 
/T 3 C 1 C 3 P23(T2<T3 of 


for (jj > 0 (j = 1,2,3), —1 < pi 2 < 1 and 1 + 2pi2/0i3/923 — P 12 ~ P 13 ~ P 23 > 0- From Kotz and 
Nadarajah (2004), the distribution of W has density 


/M = 


r((q + 3)/2) 
(a7r) 3 / 2 |E| 1 / 2 r(a/2) 


1 + 


(w — rj)'T, 1 (w — 77 ) J ( Q+3 / 2 ) 


a 


w € 


We use a cylindrical coordinate W = (X, Xi, X 2 )', where X\ = R cos© and X 2 = R sin 0 with 
R > 0 and 0 < 0 < 2ir, and consider the conditional distribution of (X, 0)' given R = r, which 
provides a distribution on the cylinder. Define new parameters as 

2 2 

p(9) = p + A cos (9 — zz), r 2 = -—^ 5 -, 

1 _ P 23 

cos p\ = r(bip 2 ~ M 3 ), «a sin^i = r(b 3 rj 3 - M 2 ), 

* n r 2 (6 3 -6i) * . r 2 b 2 

n 2 cos Zp 2 = -—-, n 2 sm ^p 2 = — 

with —00 < p < 00 , A, ref, ref > 0, 0 < 17 , pi < 27r and 0 < p 2 < 7r, where 

P = m + 017/2 + 02773, P 2 = 1 + 2/?i2/3l37'23 “ Pl2 “ Pl3 “ P23> 


Ol = 


Cl P13P23 ~ P 12 


o 2 = 


Cl P12P23 — P 13 

-t 9 ‘ 

C 3 1 - 

P23 


Then we have 

1 


C2 1 - p\ 3 

bl c 2 (l-p 2 3 )’ 62 c 2 c 3 (l - p| 3 ) ’ 

A cos i7 = —air, Asinz/ = — a 2 r. 


1 


b 3 = 


o|(i -Pi 3 y 


^(w - 77)'E 3 (u; - 77 ) = - P^)} 2 - «* cos(0 - pi) - k 2 cos 2(9 - p 2 ) + d, 

where d = ( 6 i? 7 2 + b 3 p 2 — 262 ^ 72 ^ 73 )/2 (> 0). The conditional probability density function f(x,9\r) of 
(X, 0)'|(ii = r) is represented as 


f{x,9\r) 


= C 


-1 


1 2 

1 + ^2 {x — p{0)Y - ki cos (6 — pi) — K 2 cos 2(9 - p 2 ) 


-{a+3)/2 


(1) 
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where K\ = re*/ 7 , k 2 = ^ 2 / 7 , o ’ 2 = 7 t 2 , 7 = a + d (> 0 ) and the normalizing constant is 



( 2 ) 


Here F 4 denotes Appell’s double hypergeometric function (Gradshteyn and Ryzhik, 2007, 9.180.4) 
defined by 



with Pochhammer’s symbol 



c(c+ 1 ) ■ ■ ■ (c + j ~ 1 ), J >1 

1, J = 0 


and B the beta function. The distribution with density function JT|) has nine parameters a, a > 0, 
—00 < /i < 00 , ki,K 2 ,\ > 0, 0 < v. < 27T and 0 < p 2 < tt with restriction + k 2 < 1- The 
resulting distribution should be called the generalized f-distribution on the cylinder. Relationships 
between the new and original parameters are: 

(a) 772 = 73 = 0 47 ki = 0. 

(b) K 2 = 0 47 p 23 = 0, (72 = cr 3 . 

(c) A = 0 77 P \2 = P 13 = 0. 

As being introduced in (JTJ) , the parameter a was assumed positive. However, the density is still 
valid when the parameter space of a is extended to a > — 1. Note that a corresponds to the degrees of 
freedom in the generalized f-distribution on the cylinder (JT|) and the parameter a determines the degree 
of concentration around the mode. The proposed distribution with density dH) has 9 parameters. The 
roles of the parameters are as follows: p 2 ) and p are location parameters of © and X, respectively, 

and (T is the scale parameter of X. A and v in p(6) determines the degree of correlation between X 
and 0, and A controls the skewness of X. k± and k 2 determine the morality of the distribution, 
discussed in Section 3.3. a controls the tail weight of the distribution. To see the role of a and A, 
we provide contour plots of the proposed density (H|) in Figure 1 for some combinations of a and A. 
Other parameters are specified as p = v = pi = p 2 = 0, u 2 = 1, k\ = 0.1, k 2 = 0.4. The column (I) 
in Figure 1 illustrates the interpretation of a: as a increases, the tail weight of the density increases. 
The role of A as controlling the skewness of the distribution can be seen from the column (II) in Figure 
1. The interpretation of k\,k 2 are discussed in Section 3.3. 

Under the assumption 7 = {a + 3)/2, if we use reparametrization 1/ifj = —(a + 3)/2, — k\ = 
tanh(K^*'0), — k 2 = tanh(«; 2 *' i /’)) then the density (JT|) has a similar expression to the family of distri¬ 
butions introduced by Jones and Pewsey (2005): 



i/ip 


The restriction a > — 1 in m is changed into — 1 < if) < 0 . 
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(I) 

(a) (b) (c) 





(a) 



(II) 


(b) 



(c) 



Figure 1: Contour plots of density (JTJ) for: (I) A = 0 and (a) a = —1, (b) a = 10, and (c) a = 20; (II) 
a = 6 and (a) A = 0, (b) A = 0.5, and (c) A = 1. The other parameters are set as [i = v = mi = //2 = 
0, a 2 = 1, k\ = 0.1, K2 = 0.4. 


3 Properties 


3.1 Special cases 


1. Under the assumption 7 = a/2 in (P), letting 7 (= a/ 2) —>• 00 , we have an extension of the 
distribution by Mardia and Sutton (1978). Its joint probability density function (Kato and Shimizu, 
2008) is 


f{x,6) = C 1 1 exp 


{x-KO )} 2 

2 T 2 


+ cos (0 — /il) + cos 2(9 — 1 ^ 2 ) 


(3) 


with the normalizing constant 


Ci = (2vr) 3 / 2 r 


/oK)Io(^) + 2^/ i (^)/ 2i K)cos2j(w - M2) 

l=i 


(4) 
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Here I j denotes the modified Bessel function of the first kind and order j given by 


1 /*27T 00 

Ij{z)= 2 Wo -We^d^ S 


(iH -c. 


r z£ r(r + j + l)r! 


2. When K 2 = 0, m reduces to 


f(x, 6) = C 2 1 1 + {x - /r(0 )} 2 - ki cos(0 - 7 x 1 ) 


1 —(a+3)/2 


( 5 ) 


The normalizing constant is represented as 


C '2 — 2v / 27T(jB(1/2, a/2 + 1 ) 2^1 + 1 , 1 ; re? 


using the Gauss hyper geometric function 2 -^ 1 - If we replace rei = reJ/ 7 , cr 2 = 7 r 2 and let 7 = a /2 —>• 
00 , we have the distribution proposed by Mardia and Sutton (1978) and the constant C 2 tends to 
(27 t) 3 / 2 t/o(ki)- This agrees with the fact that Ci = (27 t) 3 / 2 t/o(kJ) when re| = 0 in (|3|). 


3.2 Marginal and conditional distributions 

The marginal distribution of 0 is 

fe(0) = Cq 1 !! - kicos(0 - n 1) - K 2 cos 2(6 - ^ 2 )} _a/2_1 , (6) 


where 


■ f 1 CK o o 

C 0 = 2vr <J F 4 + 1, 1,1 ;re?,re| 


cos 2 j(fi 2 - Mi) 


, 9 (a/2 + l)3j / «i \ 2 -? /« 2 V 

( 2 j)!j! V 2 J UJ 

„ /a + 67 1 a + 67 „ „. „ _ 9 9 

x -^4 ( ^-b 2 > ^-b 1 ) 2 j + 1 , J + 1 ; Kr, «2 


(7) 


The distribution with density d 6 j) is a member of the generalized t-distributions on the circle proposed 
by Siew et al. (2008), and is possibly bimodal and asymmetric. Cosine and sine moments of the gener¬ 
alized ^-distributions are given in their paper. The generalized f-distributions include the generalized 
von Mises distribution (cf. Yfantis and Borgman, 1982) as a special case. As another special case 
when K 2 = 0 in (pQ), the marginal distribution of 0 belongs to the family of symmetric distributions 
by Jones and Pewsey (2005). Note that the marginal density dH) is independent of A, /x, <r and v which 
are the parameters of the proposed density (JT]) . On the other hand, the marginal distribution of X 
does not have a closed form in general. When A = 0, we can obtain the marginal distribution of X in 
a closed form given by 


fx(x) = C~ 1 D(x) 


1 + 2^ (I 


P-? 


— (a+3/2) 


where C is defined by © and D(x) is obtained by replacing a and «i (i = 1,2) with a + 1/2 and 
Ki/{ 1 + (x — m) 2 /(2a 2 )}, respectively, in Cq defined in ([7]). Note that this density is symmetric about 
M- 
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In Q, the conditional distribution of X given 0 = 0 has the generalized ©density function provided 


by 

f f jm = r i _ (x-/r(g)} 2 _ 

x l® \ 2cr 2 [1 — {«i cos(@ — ^i) + K 2 cos 2(0 ~ ^ 2 )}] 
with the normalizing constant 


-(a+3)/2 


C x |e = \/2aB Q, E + 1^ [1 - {ki cos(0 - Hi ) + k 2 cos 2(0 - /u 2 )}] 1/2 • 


( 8 ) 


The conditional distribution of 0 given X = x is a member of the generalized ©distributions on the 
circle given by 

/(0 \X = x) OC {1 — K\(x) cos(0 — Hi) — K 2 (x) COs2 (0 — H2)}~( a+3 ^ 2 , 


where Ki(x) = Ki/{1 + (x — /x) 2 /(2(j 2 )}, * = 1,2. As shown in Siew et al. (2008), the mean direction of 
© with density © depends on /Cj, i = 1,2. Thus the conditional mean direction E(Q\X = x ) depends 
on x through Ki(x). The result that the conditional distributions of X|0 = 0 and 0|X = x are the 
generalized ©distribution comes from the fact that the conditional distribution of a multivariate t 
distribution is again a multivariate t distribution (Joe, 2015). 


3.3 Modality 

We consider the modality of the distribution with density ©• The mode x* of (ED, whenever the value 
of 0 is specified, is x* = /x(0). Similar to Siew et al. (2008), we discuss maximization of the function 
m(0) = cos(0 — hi) + cos 2(0 — H 2 ) with respect to 0. The solution 0* of an equation 

K\ sin(0 — hi) + 2k 2 sin 2(0 — H 2 ) = 0 (9) 

is a value which maximizes m(0) if the sign of h(6*) is positive, where 

h{6) = cos(0 — Hi) + 4^2 cos 2(0 — H 2 )- 

Equation 0 can be solved numerically for any combinations of Hi an d H 2 ■ Without loss of generality, 
we let Hi = 0- Then 0 has closed form solutions when H 2 = 0,7 t/ 2,7 t /4 and 37r/4. The results 
are given in Table 1, where 0o = arccos{Ki/(4«2)}, 0i = arcsin{(— k\ + a/k^+”32k|)/(8k 2)} and 
02 = arcsin{(Ki + a/k 2 + 32ac|)/( 8/^2)}- See Yfantis and Borgman (1982) for more discussion as to the 
solutions of 0. 

Contour plots and marginal density plots of the proposed distribution are given in Figure 2. Figure 
2 shows the case when = 0.5, «2 = 0.1, A = 0, and (II) «q = 0.2, K2 = 0.3, A = 1, while the other 
parameters are set as <7 = 1, h = 0, v = 7r/3, hi = ©2 = 0, a = 6. In the case of (I) in Figure 2, the 
joint density is unimodal and the marginal densities of X and 0 are symmetric. On the other hand, in 
the case of (II) in Figure 3, the joint density is bimodal and the marginal density of X is asymmetric. 
In fact, the marginal density of X is possibly skew depending on the values of the parameters, as will 
be seen in Section 3.5. The marginal density of 0 in Figure 3 is bimodal as given in Table 1. 

In the case where a unimodal distribution is desired, we put, for example, h 2 = ©1 + vr/4 with 
restriction 2 |k 2| < to get a unimodal density function from 0. This submodel can be useful for 
modeling multimodality when a finite mixture model is easier to interpret than a bimodal distribution. 
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Table 1: Summary of the modality of the proposed distribution for some values of M 2 when m = 0. 


M2 

Condition 

Modes ( x , 0) 

0 

4^2 > «4 
4«2 < fti 

(M + A cos 0), (m — A cos z/, zr) 

(m + A cos zz, 0) 

vr /2 

4«2 > «4 
4K2 < Kl 

(M + Acos( 0 o — i^), 6 » 0 ), (m + Acos( 0 o + zz), 27 t - 0 O ) 

(M + A cos zz, 0) 

7t/4 

2 k 2 > «4 

2 k 2 < zci 

(M + A cos(0i — z/), 0i), (m — A cos(0 2 — z/), zr + 0 2 ) 

(m + A cos(@i — z/), ) 

37t/4 

2 k 2 > K4 
2 k 2 < «i 

(M — A COs (02 + Z'), 7 T — 02), (m + A cos(0l + u),2tT — 01 ) 
(M + A cos(0l + Z/), 2 - 7 T — 01 ) 


(I) 


(a) (b) (c) 





(II) 


(a) 



(b) (c) 




Figure 2: (a) Contour plot of log-density, and marginal density plots of (b) the linear random variable 
X and (c) the circular random variable 0 for the proposed distribution for: (I) = 0.5, K 2 = 0-1, 

A = 0, and (II) k\ = 0.2, K 2 = 0.3, A = 1. The other parameters are set as a = 1, m = 0, v = 7 t/ 3 , 
Mi = M2 = 0, a = 6. 
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3.4 Moments and circular-linear correlation 


We consider moments of the proposed distribution. Let a mtk = E(cosmO ) and fi m , k = E(sinkO) be 
the trigonometric moments of © under replacement a/2 with a/2 — k , which are obtainable using 
the results by Siew et al. (2008). Moments of a random vector (X, 0)' having (]Tj) are given by 

E[{X - p(Q)} 2k cosm©] = Ce tk C~ 1 2 k+1 / 2 cr 2k+1 B (k + ^-k + l\ a m , k (10) 

and 

E[{X - p(Q)} 2k sinm0] = Ce^C^+^a^B (k + ^-k + lj p m>k , (11) 

where C©^ is obtained by replacing a with a — 2k in C<~> in (|T|). 

Moreover, we derive another type of moments of a random vector ( X , 0)' having © given by 
putting «2 = 0 in ©. After some calculations, the moments for nonnegative integers k (< a/2 + 1) 
and m turn out to be 


E[{X — p(Q)} 2k cos m(0 — pi)] 

2 k a 2k B (k + i, § - k + 1) (-1 ) m+ %~ k+1 P™ /2 _ k (-(1 - K 2 )- 1 / 2 ) (1 - k 2 )-(«- 2 *+ 2)/4 


( 12 ) 


B (1/2, a/2 + 1) 2 E\ (a/4 + |, a/4 + 1,1; k 2 ) (a/2 — k - m + l) m 
where P denotes the associated Legendre function (Gradshteyn and Ryzhik, 2007, 8.711.2) defined by 


P™(z) = 


{v + l){y + 2 ) • • • (y + m) 


TV 


J (^z + z 2 — 1 cos cos imp dp. 


Next, we study the circular-linear correlation R x q between X and 0 (cf. Mardia and Jupp, 1999, 
p. 245) which is defined as 

o 2 = r ls + r lc ~ ^ cs r xs r xc 

I _ r 2 ’ 

^ 1 cs 

where r xs = Corr(A", cos 0), r xc = Corr(X, sin 0) and r cs = Corr(cos 0, sin 0) are Pearson’s correlation 
coefficients. We only consider the case when (X, 0/ has density © for simplicity because the circular- 
linear correlation of ( X , 0) with density © is not feasible to compute analytically. Note that the R x q 
in case of m can be obtained by numerical calculation. A straightforward calculation shows that 


,2 A 2 U 
' xd q + \ 2 U ’ 


(13) 


where 


U = - (1 -p 2 ) sin 2 (/ri - v) + < - (1 +p 2 ) -p\ \ cos 2 (/Ti - v), 


and pi,p 2 and q denote p m = E{cosm(Q — )},m = 1,2, and q = E[{X — /i(0)} 2 ], calculable from 

m- Note that U > 0 since 1 — p 2 > 0 and (1 + P2)/2 — p 2 = Var{cos(0 — pi)} > 0. We can 
observe that R 2 e is an increasing function of A, and R 2 g = 0 if and only if A = 0. Furthermore, letting 
7 = a/2 —>• oo, it is seen that (11311 reduces to the circular-linear correlation of the distribution proposed 
by Mardia and Sutton (1978) because © goes to the Mardia-Sutton model. This is confirmed from 
the fact that q —>• r 2 and p m —> I m (n*)/Io(K*),m = 1,2 as 7 = a/2 —>• 00 . 







3.5 Skewness of marginal distribution on the real line 

For the proposed density (fTj), we derive the skewness of the marginal density of X defined as 


7i 


E 

\X-E(X )} 3 


M 

{X-E(X )} 2 

j3/2 


Straightforward calculation shows that 

(uiA 2 + g) 3 / 2 ’ ' ’ 

where v± = Var{cos(0 — v)}, V 2 = Cov[{X — /r(0)} 2 , cos(0 ~ v)\ and V 3 = £'([cos(0 — u) — E{ cos(0 — 
z,)}] 3 ), which are calculable from (1101) and dill) . We can easily obtain from (114|) that 71 — > v;$/v'^ 2 as 
A —>• 00 and 71 = 0 when A = 0. Figure 4 shows that a graph of skewness as a function of A. We see 
from Figure 4 that the marginal distribution of X of the proposed model can be left and right skewed 
according to the values of the parameters. 


(a) 



x 

Figure 3: Skewness of the marginal distribution of X 
0.5, K2 = 0, Hi = H2 = 0, a = 9 with v = it (solid), 
fj = 3, // = 0, v = 0, K\ = 0.3, Hi = 1^2 = 0, a = 
K 2 = 0.2 (dotted). 


(b) 



0 10 20 30 40 50 

X 

5 a function of A for (a) a = 3, /i = 0, k\ = 
= 7t/2 (dashed) and v = 0 (dotted), and (b) 
9 with K 2 = 0 (solid), K 2 = 0.1 (dashed) and 


3.6 Regression 

We can derive a circular-linear regression model from the conditional distribution with density (|5|). In 
fact, the conditional mean of X given 0 = 9 is 

E(X |0 = 6 ) = /r(0) = h + A cos(0 — v) 


9 















and the conditional variance of X given 0 = 9 is 


Var(X|0 = 9) 


2a 2 


a 


{1 — Kl cos (9 — /il) 


k,2 cos 2(6 — ^ 2 )} • 


Note that the conditional variance is dependent on 6 , i.e. the regression model is possibly heteroge¬ 
neous. As a reduced model, if we let 7 = a —>■ 00 , we have Var(X|0 = 6) = t 2 , which is independent 
of 6. Moreover if k\ = 0 and «2 = 0 in |S]), we have Var(A|0 = 9) = 2a 2 /a and, in this case, we 
obtain a regression model 


Xi = p + A cos (6i — v) + £i, i = 1,..., n, 

with random errors which are independent and identically distributed according to the generalized 
t-distribution. 


3.7 Maximizing /3-entropy 

The related distribution proposed by Mardia and Sutton (1978) and Kato and Shimizu (2008) can 
be characterized as the maximum entropy distribution under certain moment conditions. Also a 
maximum entropy distribution under certain moment conditions relates to the proposed distribution 
with density P). We consider the /3-entropy (see Eguchi, 2009, Section 13.2.4) defined as 

E (f) = {J f +1 (x,0)dxd6-p-l\ f 1 3(/3 + l). 

Then the maximum entropy distribution subject to constraints on the moments 

E(X 2 ), E(X cos 0), E(X sin 0), F(cosp0), E(sinpQ), p= 1,2, 


is the distribution with density 


f(x,6) oc 1 + /1 


{x-p(0)Y 

2 T 2 


+ cos (9 — p\) + K 2 cos 2(9 — P 2 ) 


VP 


(15) 


If we take /3 = — 2/(a + 3), (fl5l) gives a density related to (P). 


3.8 Parameter estimation 


We provide a method for calculating the maximum likelihood estimates of the parameters in the 
generalized t-distribution with density P). When we observe ( Xi,9i ), i = 1 ,... ,n, the log-likelihood 
function is given by 


1 Q( \ Tl 

L(ip) = -nlogB (-,- + !)- n log Cq - ^°&2 - nlogcr 


1 


x (« + 3) lo g 


i= 1 


1 2 

1 + 7 —o i x i — P - A COS (9i — v)\ — Ki cos (9i - pi) - K2 cos 2 (9i — p 2 ) 

2 ( 7 - 


where ip = (a, p, A, v, kj, pi, K 2 , P 2 , &)'■ For obtaining the maximizer of L(ip), we propose the con¬ 
ditional maximization algorithm. We first divide the parameter ip into ip = (a, ip^, ip?)', where 
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•01 = (fJ,,X,u)' and 0 2 = (^ 1 , fi\, « 2 > M 2 , Given the value of <7 and 0 2 , maximizing L(0) is 
equivalent to maximizing 


-E lo § 

2 — 1 


A + ^2 {Xj - (J, - A cos(0* 



with respect to r/q, where A = 1 — fiq cos(0* — mi) — K 2 cos 2(0* — ^ 2 )- Let x = (aq,..., x n ) r , T = 
(A ..., t' n Y for £* = (1, cos 0*, sin 0*)', (3 = (//, A cos z/, A sin z/)' and VF = diag(u;i,..., uz n ) for 


a 2 + {a;* — h — A cos(0* — zz)} 2 /2 


Using the theory of weighted regression (see Andrews, 1974), the maximizer of 0 1 given a and 0 2 can 
be obtained as 

3 = (T , WT)~ 1 T , Wx, (16) 

which deduces the maximizer r/q. Since W depends on t/q, we calculate W based on the current 
values in each iteration. 

Given r/q and 0 2 , maximizing L(0) with respect to cr 2 is equivalent to solving the following 
equation: 

( , qn-i _ {xi-n~ Acos(0*-z^)} 2 

72 (CK + o) — / - 9 - -5 

{%i ~ H ~ A cos(0* - r')} 2 + Cjcr 2 

which deduces the maximizer a. 

Finally for maximizing 0 2 under given a and r/q, we maximize 



— n log B 



n log Cq - - (a + 3) ^ log {A - «i cos(0* - /xi) - k 2 cos 2(0* - M 2 )} (18) 

i=l 


with respect to 0 2 , where A = 1 + (2cr 2 ) -1 {x* — m — A cos(0* — z/)} 2 . This maximization problem is 
quite similar to obtaining the maximum likelihood estimates of the generalized f-distribution on the 
circle (Siew et al., 2008), so that we can obtain the maximizer 0 2 given a and r/q. Note that we 
carried out the maximization of (1181) with use of numerical integration for getting the value of Cq. 

Therefore, the proposed estimation method is described in the following. 


Estimation Algorithm 

1. Determine the initial values r/A\ Set k = 0. 

2. Compute W based on and then obtain 'ip[ k+1 ' > based on (fl6l) . 

3. Compute A fc+1 ) by solving (fT71) with r/q = 0^ A+1 ' and r/> 2 = 0 2 ^. 

4. Compute 0 2 fc+1 ^ by maximizing (flHl) with r/q = 'ip[ k+1 ' 1 and a = a^ k+1 \ 

5. Set k = k + 1 and go to Step 2 (until numerical convergence). 
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4 Empirical application 


For an illustrative example, we consider a cylindrical dataset given in Johnson and Wehrly (1977) on 
the wind direction and ozone level taken at 6:00 pm at four-day intervals between April 18th and June 
29th, 1975 at a weather station in Milwaukee with 19 samples. We fitted the proposed generalized 
f-distribution and its submodels. For submodels of dU), we consider two cases, namely, ^2 = 0 (GT- 
subl) given in ([5]) and /i 2 = H i + tt/ 4 with 2\^\ < «l (GT-sub2) discussed in the end of Section 3.3. 
When we carry on the estimation algorithm given in Section 3.8, we repeat the algorithm until the 
difference between update and current values are smaller than 0.001. For comparison, we also fitted 
the member of the exponential family given by Kato and Shimizu (2008) with density 


f(x,6) cx exp ^ - 


(x-n{0)) 2 


2cr 2 


+ Ki COS (9 — fi\) + K -2 COS 2(0 — H2 ) 


which is reparametrized form of ([3]) with k* = k\, = K 2 and r 2 = cr 2 . Table [2] provides the 

maximum likelihood estimates of the parameters, AIC values, and the multivariate Kolmogorov- 
Smirnov statistics of goodness of fit testing given by Justel et al. (1997). In Table 1 of Justel et 
al. (1997), the bivariate Kolmogorov-Smirnov statistics distributions are reported. From the table, 
the upper 5%, 10% and 25% percentiles are 0.362, 0.335 and 0.292, respectively, for 20 sample cases. 
Thus the percentiles in 19 sample cases are smaller than these values, so that the fitted four models 
are not rejected with 5% significance level. Judging from AIC, we see that the two submodels of the 
proposed distribution, GT-subl and GT-sub2, gives better fits than the Kato and Shimizu distribution. 
Figure [3] shows scatter plots of the data and contour plots of the fitted densities of two submodels. 


Table 2: Maximum likelihood estimates of the parameters, the maximum log-likelihood (MLL), and 
AIC values of model JT]) (GT) and its submodels and the Kato and Shimizu model d3|) (KS) fitted to 
the data from Johnson and Wehrly (1977). 


Model 

h 

A 

V 

a 

«1 

/A 

H2 

M2 

a 

AIC 

g.o.f 

GT 

41.38 

31.14 

1.25 

68.90 

0.11 

6.28 

0.03 

1.35 

23.84 

417.79 

0.180 

GT-subl 

41.37 

31.17 

1.24 

74.81 

0.08 

0.30 

— 

— 

28.28 

240.20 

0.314 

GT-sub2 

41.01 

32.70 

1.41 

6.08 

0.49 

0.19 

0.15 

— 

-1.00 

203.34 

0.359 

KS 

41.24 

31.38 

1.20 

19.77 

1.41 

0.19 

0.35 

1.45 

— 

241.42 

0.282 


5 Conclusions 

In this paper, we derived the distributions on the cylinder based on a trivariate t-distribution. The 
derived distribution is considered as a cylindrical extension of the generalized t-distribution on the 
circle proposed by Siew et al. (2008) and includes the exponential family given by Kato and Shimizu 
(2008). We investigated some properties of the proposed distribution including the marginal and 
conditional distributions, circular-linear correlation and the algorithm for parameter estimation. We 
applied our proposed distribution to the data set of the wind direction and ozone level given in 
Johnson and Wehrly (1977), and we confirmed that the proposed distribution gave a better fit than 
the distribution given by Kato and Shimizu (2008). 
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Figure 4: Sample plot and contour plot of the fitted density function. 
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